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Abstract. The paper analyses the decay of any zero modes that 
might exist for a massless Dirac operator H := a-{l/i)'V+Q, where 
<5 is 4 X 4-matrix- valued and of order 0(|x|^^) at infinity. The 
approach is based on inversion with respect to the unit sphere in 
and establishing embedding theorems for Dirac-Sobolev spaces 
of spinors / which are such that / and Hf lie in (^^(R"^)) , 1 < 
p < oo. 



1. Introduction 

The mathematical interpretation of the stabihty of matter problem 
concerns the question of whether the energy of a system of particles is 
bounded from below (stability of the first kind) and by a constant mul- 
tiple of the number of particles (stability of the second kind). Dyson 
and Lenard made the initial breakthrough in 1967 for a non-relativistic 
model, and since then the problem has attracted a lot of attention, var- 
ious relativistic models having been intensively studied in recent years. 
As was demonstrated by Frohlich, Lieb and Loss in [3], to establish sta- 
bility, it is of crucial importance to know if the kinetic energy operator 
has zero modes, i.e. eigenvectors corresponding to an eigenvalue at 0; 
the possibility that zero modes can exist was established at about the 
same time by Loss and Yau in [3]. Subsequently, Balinsky and Evans 
showed in [I] that zero modes are rare for these problems. 

The first objective of the research reported here was to confirm that 
the results in [1] can be extended to Dirac-type operators with matrix- 
valued potentials. This then set the scene for the main goal which was 
to determine the decay rates of zero modes whenever they occur. The 
Dirac operator considered is of the form, 

H = a-p + Q, p = -iV (LI) 

where a is the triple of Dirac matrices and Q is a 4 x 4 matrix- valued 
function. Assuming that ||(5(-)||c4 € //^(M^), where || ■ He* denotes any 
matrix norm on C^, it was shown that Q is a small perturbation of 
a ■ p and hence fll.lj) defines a self-adjoint operator H as an operator 
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sum with domain (//^'^(M'^)) , the space of 4-component spinors in 
(L^(]R^))*^ with weak first derivatives in (L^(M^))^. The technique in 
P readily apphed to (11.11) to meet the first objective and yield the 
following result: 

Theorem 1. Let ||Q(-)||c4 e L^{M.^). Then Ht := a ■ p + tQ,t e 
can have a zero mode for only a countable set of values oft. Moreover 

nul(H) < const. / |Q(x)|Mx, 

where nul(H) denotes the nullity of H , i.e. the dimension of the kernel 
ofH. 

It is proved in [B] (see also [7j) that if ||Q(x)||c4 = 0(< x >"''), p > 1, 
where < x >:= -^/l + jx^, then any zero mode is 0(|x|^^) at oo and 
there are no resonances (defined as a solution say which is such that 
< X ip e (L^(]R^))^ for some s > 0) if p > 3/2 : note that such 
a Q satisfies the assumption of Theorem 1. In view of this result we 
concentrate our attention on more singular potentials Q, namely ones 
which satisfy 

||g(x)||c4 = 0(|x|-i), XG5J, (1.2) 

where Bi is the unit ball, centre the origin in M^, and denotes 
its complement. We actually consider weak solutions of Hip = in 
(L^(i?5^))^ . Hence the behaviour of Q at the origin, and indeed within 
Bi, does not feature. Our main result is given in Theorem 4 below. Our 
approach is very different to that in [6j and is based on two techniques. 
In the first we use inversion with respect to Bi to replace the problem 
in BI by an analogous one in Bi. The second part involves establishing 
Sobolev-type embedding theorems for the spaces EI^''f'(f2) defined as the 
completion of [C^(i7)]^ with respect to the norm 

\\f\Un:=y^i\{a-p)f\'+\f\nd^y' . (1.3) 

For p = 2 these are just vector versions of the standard Sobolev spaces, 
but we need cases p 7^ 2 for which we were unable to find any ap- 
propriate results in the literature. The results we give in section 3 
use a method of Ledoux in ^ to derive a weak inequality and then 
embedding properties of Lorentz spaces on bounded domains. 

2. Reduction by inversion 

We recall that in (11.11) . a = (cti, 02, eta) are Hermitian 4x4 matrices 
satisfying 

ajak + akOj = 2Sjkh, j, k = 1,2, 3, 
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where /4 is the unit 4x4 matrix, and we choose 

02 ff, 



"'=1^; oi )• (2-1) 

where the aj are the Pauh matrices 

-^(°o')' 

and O2 is the 2x2 zero matrix. From (11.21) we have that the 4x4 
matrix- valued function Q has components qjk,j,k = I,-- - ,4 which 
satisfy 

|g,-fc(x)| <C|x|-\ |x|>l. 

For C^-valued functions /, g, measurable on an open set C and 
z,w & C^, we shall use the following notation: 

4 

< z,w > := ''^^ZjWj, \z\:=<z,z>^ 



{f,9)n ■= < /W,^(x) > rfx 
Jn 



i/p 

i/(x)rrfx' 



Thus {f,g)n is the standard inner-product on (L (fi)) and ||j ||!^ 
||/||2n the standard norm: when f2 = we shall simply write {f,g) 
and i|/||. 

Let ip he a weak solution of HiIj = in \ where Bi is the 
open unit ball centre the origin and Bi is its closure: hence for all 
<jye{C^{R'\B,))' 

I := j < if7/'(x),0(x) > cix = 0. (2.3) 

Inversion with respect to Bi is the involution Inv : x i-^ y, y = x/|xp, 
and for any function defined on \ the map M : 1^ := 
0o Inv~^ is such that 0(y) = 0(x) and yields a function on Bi. Hence 

e {C^{M?\Bi)f means that G {C^{Bi \ {0}))^ The inversion 
gives 

M{(a-p)V'}(y) = |yP(/5-p)^(y), (2.4) 
where [3 = {^1,^2,^3) and 

3 
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where Sj is the Kronecker delta function. It is readily verified that the 
matrices (3k{y) are Hermitian and satisfy 

Also there exists a unitary matrix X{y) such that X G C°^(]R'^ \ {0}) 
and for all y 7^ 0, 

X{yr'Pk{y)Xiy) = k = l,2, 3. (2.5) 

Setting LJ := y/|y| it is easy to verify that these conditions are satisfied 
by 



where 



—UJ2 + iuJi —iuj-s 



Let ipiy) = '~^(y)^(y)- Then from (12.41) we have 

M{(a ■ p)^}(y) = \y\'Xiy) {{a ■ p)^{y) + r(y)^(y)} , (2.7) 
where 

(2.8) 



Y{y) = j2c^kXiy)-' f-^T^X(y; 

k=i \ 



Also, a calculation gives that the Jacobian of the inversion gives dx = 
\y\''dy. 

Returning now to (12. 3p . and with 0(y) = — X(y)$(y) , the inversion 
yields 

/ = / <|y|2X(y){(«.p)vI/(y)+F(y)vI/(y) 

JBi 

- g(y)X(y)M/(y)} , X(y)<|.(y) > \y\-'dy = 
for all $ e {C^{Bi \ {0}))% which can be written as 

/ = / < (« . p)vl/(y) + Z(y)xl/(y), |yr^$(y) > dy = (2.9) 

JBi 

where 

Z{y) = Y{y) - \yr'X{y)-'Q{y)X{y). (2.10) 
Equivalently, we can remove the factor |y|"Mn to give 

/= / <(«-p)^(y) + Z(y)^(y),$(y)>rfy = (2.11) 

JBi 

for all $ e {C^{Bi\{0}))\ 
Let C e C°°(M+) satisfy 



for < t < 1 

1 for t > 2 
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and for y G set Cn(y) = C('^lyl)- Then 

VCn(y)=<(n|y|)^ 

and so 

|vay)| = o(nx[i,i](|y|)^ 

where xi denotes the characteristic function of the interval /. We then 
have from dm]) , now for all $ G {C^{Bi)f , 

< (a ■ p)M/(y) + Z(y)^(y), Cn(y)$(y) > dy 

Bi 



< Cn(y) {(« ■ P)^(y) + ^(y)^(y)} , $(y) > dy 

Bi 

= 0. (2.12) 

In J, 

C„(y)(a ■ p)M/(y) = {a ■ p)(Cn^)(y) - [(« • p)Cn](^)(y) 

and 

K(y) := [{a ■ p)Cn] = 0(nxri 2,). 
Therefore, as n ^ cxd, 

|/ <Vn{y)^y),^y)>dy\ < uW^Wl^^b^) X[i.i]dy 

= 0(n-2)^o. 

Also 

f <[(a-p)+Z(y)]C„(y)vi>(y),$(y)>f/y 

J Bi 

= [ < Cn(y)^(y), [(« ■ P) + ^(y)]<^'(y) > ^^y 

JBi 

^ / <vI/(y),[(a.p) + Z(y)]$(y)>rfy. 

JBi 

We have therefore proved that 

/ <[(a-p) + Z(y)]vI/(y),$(y)>rfy = (2.13) 

JBi 

for all $ G (C^(i?i))'^ . In other words 

[(a ■ p) + Z(y)]vl/(y) = (2.14) 
in the weak sense. From fl2.10p it follows that 

Il^(y)||c4<qyr^ (2.15) 



6 A. BALINSKY, W. D. EVANS, AND Y. SAITO 

3. DiRAC-SOBOLEV INEQUALITIES 

Let If'P{n),l < p < oo, denote the completion of [C^{n)]^ with 
respect to the norm 

i,p;o:={^(|(«-p)/r+|/ncix| ^\ (3.1) 

We also use the notation 

©:=(«■ p)2, Ft : e-*, t > 0. (3.2) 

Then, 

(D/), = -A/,, (P,/)^. = e-*^/,, J = 1,2,3,4, 
where {e~*^}t>o is the heat semigroup. Furthermore, for all t > 0,x G 

and so 

Note that if 7^ M" we put any / G ]H[^'P(f^) to be zero outside f2 and 
hence is in H^'^ = M^'PfM"). Define 



iB^in) :=sup{r"/2|pj|^.^}, (3.5) 

t>o 

where |Pt/|oo;n := sup^gn |Pt/(x)|, and denote by B"'{Q) the comple- 
tion of [C^(f2)]'^ with respect to || ■ ||BQ(f^). 

Our main theorem in this section introduces the weak-L'^ space on 
fl, written, L'^'°°{Q), which is defined by 

\q,oo;n ■ sup {u'^X{\f\ > u)} 



u>0 

where A denotes Lebesgue measure and A(|/| > u) stands for the mea- 
sure of the set in Q on which |/(x)| > u. 

Theorem 2. Let 1 < p < q < 00 and let f be such that || (a • p)/||p ^ < 
00 and f G B^^^^~^\il), for 9 = p/q. Then we have for some constant 
C > 0, 

\\f\Uoo;n<C\\{a-p)f\\l (3.6) 

Proof. For simplicity of notation, we suppress the dependence of the 
norms and spaces on Q throughout the proof. It is sufficient to prove 
the result for / G {C^{Q))'^ . The proof is inspired by that of Ledoux 
in [3]. By homogeneity, we may assume that < 1, and so 

for all t > 0. Thus, on choosing t„ = u'^(^~^^^^ it follows that 

Moo < u. (3.7) 
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This gives that |/| > 2u imphes that \ f - FtJ'l > \ f\ - \Ftf \ > u and 
consequently 

u'^K\f\>2u) < n^A(|/-nj >n) 

< W^^^ j |/-Pi„rrfx. (3.8) 

Since 

^P,/ = (« ■ p)¥,/, Po/ = /, 

in view of the analogous result for e~*^ on each component of /, it 
follows that ^ 

Jo 

and hence for all g G [C^{Bi)]'^ , 
<gJ-FJ>d7c = 



(^j ^<{a- p)P,^, {a ■ p)/ > rfx^ ds 



< ||(a-p)/||p / \\{a-p)¥sg\\p:ds, 

Jo 

where p' = p/{p — ^) for p > 1 and p' = oo otherwise. From (13.51) we 
have 

1 ^ /* d ^ 

= (i^^l^'^^t"-^"-^)]^^^)'^^- 

On using Young's inequality for convolutions, this yields 

1 If 

||(a-p)P,^||p/ < (4^g-)3/2 ^ Ls^^'^^^ ''^^ 

< Cs-i 11^711,, (3.9) 
for all p G [1, oo). We therefore have 

I J <gJ-¥J>d^\<Ct^\\{a-p)fmg\\p' 

and thus 

\\f-FJ\\,<Ct^\ia-p)f\\,. 
On substituting this in (13.81) we have 

u'^X{\f\>2u) < Cu'^~nPjma-p)fr^ 
= C\\ia-p)f\\l 
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since q — p + p{6 — l)/6 = 0, whence the result. 



□ 



In the following corollary, the notation indicates that integration is 
over Bi. 

Corollary 1. Let 1 < p < q < oo,r := 3(^ — 1) G and f G 

EI^'^(i?i). Then we have that for any k G (0, g) and 9 = p/q, there 
exists a positive constant C such that 

1-6 



<cii(«.p)/ii^,^ji/ii:- (3.10) 

Proof. All norms in the proof are over Bi. From (13. 5p . for any r G [1, oo] 
and with r' = r/(r — l),r>l, 

(47rt)2 \Jb, ) 



note that this holds also for r = 1. Hence, 

= C\\f\\,. 

if 2{e-\) ^ 'h ~ which is true if r = 3(g/p — 1) since Q = p/q. 
Since r < p and / G U'{Bi) it follows that / G U{Bi) and hence 
/ G From fl3.6l) we therefore have 

||/||,,oo<C||(a-p)/||J||/||^^). (3.11) 

But for Lorentz spaces L^'** on a set Vt of finite measure, we have the 
continuous embeddings (see [2], Proposition 3.4.4) 

L'i^'iSl) ^ L'^'™ (3.12) 

if < A; < g < oo, < s, m < oo. In particular, with s = oo,m = k, 
and recalling that L^''^ = L'^, we have for < A; < g < oo, 

and so 

ll/IU<C||/ll,,oo. 

The corollary follows from (13.111) . □ 

Corollary 2. Let p G [l,oo),A; G [l,p(p + 3)/3) anc? / G HI'P(Si). 
r/ien i/iere exists a positive constant C such that 

<C||(«-p)/|Ub,. (3.13) 
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Proof. For k G we choose in Corollary 1, q = p{k + 3)/3 and 

r = k to deduce fl3.13l) from fl3.10p . When k G (p, 4p/3) we choose any 
q G [4p/3,p(p + 3)/3], so that q > k > p and r = 3{q/p — 1) G [l,p]- 
When k G [4p/3,p(p + 3)/3) we choose any q G {k,p{p + 3)/3] so that 
q > k > p and r G In both the last two cases we also have r < k 

and hence 

\\f\\r,B.<C\\f\\,,n.. 
This yields (l37[3D from (l310D . □ 

4. Estimate for zero modes 
Let ip be such that ip, (a ■ p)^ G L^{B^), := \ Bi, and 

(a-p)^(x) = -Q(x)^(x), (4.1) 

where 

||Q(x)||c4 = 0(1x1-1), 

and 

i^e{L\Bt))\ _ 

If G Ci(]R^) is 1 in the neighbourhood of oo and is supported in M^\i?i 
then Oip has similar properties to those above. Hence we may assume, 
without loss of generality, that ip is supported in R'^ \ Moreover, 

/ |xn(a ■p)V^(x)prfx < cx). (4.2) 

On applying the inversion described in section 1, and using the notation 

(M^)(y) =:V^(y) =: -X(y)vl/(y) 
where |X(y)| x 1, we have from ip G L^(B^) that 

/ l^(y)l'ni<^ (4-3) 
Jbi \y\ 

We also have from (12. 7p that 

(M(a ■ p)^) (y) = |y|2 {X(y) [{a ■ p)vl/(y) + Y{y)^iy)]} , (4.4) 
where V(y) is given in (12.81) and is readily seen to satisfy 

r(y)llc^^i/|y|- 

Let ^(y) = |yp$(y). Then 

(a-p)vI/(y) = |y|2|(a.p)<|.(y) + 
and so from (14.41) 

\M{a ■ p)^(y)| X |yr |(« ■ p)<I.(y) + O (^|^) } • (4-5) 
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Hence from (14. ip 









Mi 







(a-p)$(y)+0 



fmy)\\ 



Since, from (4.3), 



V |y| J 

I |$(y)p^<oo 



2 



^ < oo. (4.6) 



it follows from (14. 6 p that 

/ |(a-p)$(y)pc/y < oo. 

J Ex 

Hence, $ G M^^'^{Bi) : recall that we may assume that ip is supported 
in \ Bi and hence $ is supported in Bi. By Corollary 2, we have 
that $ e L*^(5i) for any k e [1, 10/3). 
We therefore have the preliminary result 

Theorem 3. Let ip ^ L'^i^f) he a solution of (^T7]) with Q{x) = 
0{\x.\^^) in B^. Then, for any k G [1, 10/3), 

v;(x) = ixrv(x), 

where 

|(/)(x)|'=|xr^rfx< oo. 



However, we can repeat the analysis that led to Theorem 3 to improve 
this result. From '^{y) = |yp$(y) and (2.14) it follows that 

(a-p)$(y) = -Z«(y)$(y), (4.7) 

where 

Z« = Z(y)-2^|y|-2(«-y) = 0(^). 

|y| 

Let $(y) = |y|*$^^''(y). Then we have 

(a-p)<l>(2)(y) = -Z(2)(y)$(2)(y), 

where 

By Holder's inequality, for 1 < p < A;, /c G [1, 10/3), 
f |z(2)(y)$(2)(y)|^rfy < c [ \\y\-^'^'^Hy)\'dy 

JBi J Bi 

< / mytdy] / \yr^'^^)Wik-,)^\ 



Bi / \JBi 



< oo 

if 
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Choose p = 1 and let < t < 11/10. Then fl4.8p is satisfied for some 
k e [1, 10/3), depending on t. We therefore have that G EI^'^(i?i) 
if < t < 11/10, and hence, by Corollary 2, that <l>(^) G L''{Bi) for any 
s G [1,4/3). The following result has consequently been proved. 

Theorem 4. Suppose that'ip G L?'{Bl) satisfies {(a-p) + (5(x)}^/'(x) = 
in Bf. Then for any < t < 11/10, = |x|^^^*(/)(x), where 
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